Astronomy & Astrophysics manuscript no. n 1377 'aalto 


©ESQ 2013 


March 8, 2013 





O 



Winds of change - a molecular outflow in NGC 1377? 
The anatomy of an extreme FIR-excess galaxy 

S. Aalto 1 , S. Muller 1 , K. Sakamoto 2 , J. S. Gallagher 3 , S. Martin 4 , and F. Costagliola 1 

1 Department of Earth and Space Sciences, Chalmers University of Technology, Onsala Observatory, SE-439 94 Onsala, Sweden 
e-mail: saalto@chalmers . se 

2 Institute of Astronomy and Astrophysics, Academia Sinica, P.O. Box 23-141, Taipei 10617, Taiwan 

3 Department of Astronomy, University of Wisconsin-Madison, 5534 Sterling, 475 North Charter Street, Madison WI 53706, USA 

4 European Southern Observatory, Alonso de Cordova 3107, Vitacura, Casilla, 19001, Santiago 19, Chile 

Received x; accepted y 



a,: 



£ ■ ABSTRACT 

Aims. Our goal was to investigate the molecular gas distribution and kinematics in the extreme FIR-excess galaxy NGC 1377 and to 
address the nature and evolutionary status of the buried source. 

Methods. We use high (0."65 x 0/'52, (65 x 52 pc)) and low (4/'88 x 2."93) resolution SubMillimeter Array (SMA) observations to 
image the 12 CO and 13 CO 2-1 line emission. 

Results. We find bright, complex 12 CO 2-1 line emission in the inner 400 pc of NGC 1377. The 12 CO 2-1 line has wings that are 
tracing a kinematical component which appears perpendicular to the component traced by the line core. Together with an intriguing 
' X-shape of the integrated intensity and dispersion maps, this suggests that the molecular emission of NGC 1377 consists of a disk- 

outflow system. Lower limits to the molecular mass and outflow rate are M out (H 2 )>l x 10 7 M© and M>8 M Q yr _1 . The age of the 
proposed outflow is estimated to 1.4 Myrs, the extent to 200 pc and the outflow speed to V out =140 km s _1 . The total molecular mass in 
the SMA map is estimated to M tot (H 2 )=1.5 X 10 8 M© (on a scale of 400 pc) while in the inner r-29 pc the molecular mass is M core (H 2 ) 
= 1.7 X 10 7 M Q with a corresponding H 2 column density of Af(H 2 )=3.4 x 10 23 cm -2 and an average 12 CO 2-1 brightness temperature 
of 19 K. 13 CO 2-1 emission is found at a factor 10 fainter than 12 CO in the low resolution map while C 18 2-1 remains undetected. 
. We find weak 1 mm continuum emission of 2.4 mJy with spatial extent less than 400 pc. 

Conclusions. Observing the molecular properties of the FIR-excess galaxy NGC 1377 allows us to probe the early stages of nuclear 
activity and the onset of feedback in active galaxies. The age of the outflow supports the notion that the current nuclear activity is 
young - a few Myrs. The outflow may be powered by radiation pressure from a compact, dust enshrouded nucleus, but other driving 
mechanisms are possible. The buried source may be an AGN or an extremely young (1 Myr) compact starburst. Limitations on size 
and mass lead us to favour the AGN scenario, but further studies are required to settle the issue. In either case, the wind with its 
implied mass outflow rate will quench the nuclear power source within a very short time of 5-25 Myrs. It is however possible that the 
gas is unable to escape the galaxy and may eventually fall back onto NGC 1377 again. 

*o ' 
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^vq . 1 . Introduction well as insights into how the infrared and radio emission is reg- 
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, , ulated in galaxies. 

A small subset of galaxies deviate strongly from the 

well-known radio-to-FIR correlation through having ex- 
cess FIR emission as compared to t he radio (a > 3; 1 .1 . The extreme FIR-excess galaxy NGC 1377 



cess f ik emission as comparea to tn e radio (g > i; 
4=log[FIR/3.75xl0 12 Hz]/S„(1.4GHz) (IHelou et alJ Il985lf ). 
^ I There are several potential interpretations ofthe excess including NGC 1377 is a member of the Eridanus galaxy group at an esti- 
very young synchrotron-deficient starbursts or dust enshrouded mated distance of 21 Mpc (1"=102 pc) and has a far-infrared 
Active Galactic Nuclei (AGN). AGNs obscured by nuclear dust luminosity of L mR = 4.3 x 10 L ( |Roussel et al. | |2003j). In 
are likely also in the early stages of their evolution where nu- stellar light, NGC 1377 has the appearance of a regular lenticu- 
clear material has not yet been dispersed and/or consumed to lar § alax y ( |de Vaucouleurs et al.| | l99l|) , with a diameter of 1'.8, 
feed the growth of the black hole. The FIR-excess galaxies are a l^ge scale inclination of 60° and a major axis position an- 
rare - iRoussel et all d20f)l find that they represent a small frac- gle (PA) of 92° (asdenved from the K-band image). However, 
tion ( 1 %) of an infrared flux-limited sample in the local universe, i Heisler & Vadei] <|1994j) reported the presence of a faint dust lane 
such as the IRAS Faint Galaxy Sample. Despite their scarcity, a l° n § the southern part of the minor axis, perturbing an oth- 
which likely is an effect ofthe short time spent in the FIR-excess erwise featureless morphology. NGC 1377 is a Sixty Micron 
phase, these objects deserve careful study. Their implied youth Peaker ( SMp ) meaning that its IR SED peaks near 60/rni. SMPs 
provides an ideal setting to better understand the initial condi- often found t0 have peculiar morphologies and are either 
tions and early evolution of starburst and/or AGN activity, as classified as H II region like or Sy 2 - where the Sy 2 are some- 
wha t more common (e.g. IHeisler & Vad er 1994; lLaureijs et al.l 
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NGC 1377 is the most extreme example known so far of an 
FIR-excess galaxy with radio synchrotron emission being defi- 
cien t by at least a factor of 37 with respect to normal galax- 
ies dRoussel et al.1 12003L l2QQ6h (corresponding to q > 3.92). 
Interestingly, H II regions were not detected through near- 
infrare d hydrogen recombina tion lines or thermal radio contin- 
uum (Roussel et al. 2003, 2006). The presence of cold molecular 
gas in the center of NGC 1377 is evident through bright single- 
dish 12 CO 1-0 and 2-1 line emission detected by Rou ssel et al.l 
(2003). Deep mid-infrared silicate absorption features suggest 
that the nucleus is enshrouded by large masses of dust (e.g. 
ISpoon et aD l2007) potentially absorbing all of the ionizing pho- 
tons. The extremely high obscuration aggravates the determina- 
tion of the nature of the nuclear activity. 

iRoussel et all (l2006h propose that NGC 1377 is a nascent 
(ft 1 Myr) opaque starburst - the radio synchrotron deficiency 
caused by the extreme youth (pre- supernova stage) of the star- 
burst activity where the y oung star s are still embedded in their 
birth-clouds. In contrast, Imanishi| (120061) argue, based on the 
small 3.3 jim PAH equivalent widths and very red L-band con- 
tinuum, that NGC 13 77 harbours a buried AGN. Furthermore, 
llmanishi et all d2009h find an HCN/HCO + /=l-0 line ratio ex- 
ceeding unity, which they suggest is evidence of an X-ray 
Dominated Region (XDR) surrounding an AGN. They explain 
the lack of radio continuum through suggesting the presence of 
a large column of intervening material causing free-free absorp- 
tion. Note, however, that the HCN/HCO + line ratio determina- 
tion of the XDR contribution is uncertain a nd is based on small 
differences in the ratios (see e.g. Fig. 5 in iKrips et aD ([2008)). 
Furthermore, the lack of a thermal radio continuum detection 
also put limits on the possibility of nuclear free-free absorption. 

To further our understanding of the nature of NGC 1377 
we obtained high resolution 12 CO and 13 CO 2-1 data with the 
SubMillimeter Array (SMA) in Hawaii. As NGC 1377 remains 
undetected in HI and Ha a ccording to reports in the literatureQ 
dRoussel et al.l 120031 120061) the molecular lines seem to be 
an important way to study the gas properties and dynamics. 
Our goal was to search for clues to the extreme FIR-excess 
of NGC 1377 through the distribution and kinematics of the 
molecular gas. We found evidence of a molecular disk-outflow 
system and we we use its properties to address the nature of the 
buried source. 

In Sect. [2] we present the observations and results which are 
discussed in terms of the structure and power of a molecular 
disk-outflow system in Sect. 13. II and Sect. 13.21 In Sect|33] pos- 
sible power sources of the outflow are presented and in Sect. 13.41 
we briefly compare with molecular outflows in the l itera ture. The 
properties of the nuclear gas are discussed in Sect. 13.5 l and evo- 
lutionary implications for NGC 1377 are discussed in Sect. 13.61 

2. Observations and results 

NGC 1377 was observed with the SubMillimeter Array (SMA) 
on 2009 July 17th, in the very-extended configuration (8 anten- 
nas), and on 2009 October 12th, in the compact configuration 
(7 antennas). The phase center was set at #=03:36:39.10 and 
£=-20:54:08.0 (J2000). During both nights, the zenith atmo- 
spheric opacity at 225 GHz was near 0.1, resulting in system 

1 Gallagher et al (2012 in preparation) obtained a spectrum with the 
Southern African Large Telescope (SALT) that shows emission for Ha 
and [N II], and we therefore assume that these lines, along with the 
previously measured [S II] double, are weak but present in NGC 1377. 



temperatures in the range 100 - 200 K depending on source ele- 
vation. 

The heterodyne SIS receivers were tuned to the frequency 
of the 12 CO 2-1 transition at 230.538 GHz in the upper side- 
band, while the 13 CO 2-1 transition was observed in the lower 
sideband. The correlator was configured to provide a spectral 
resolution of 0.8125 MHz. 

The bandpass of the individual antennas was derived from 
the bright quasars J1924 - 292 in July, and 3C454.3 in October. 
The primary flux calibration was set on Callisto and Uranus, 
respectively. The close-by quasars J0423-013 (~4 Jy at 1 mm) 
and J0334-401 (-1.3 Jy) were observed regularly for complex 
gain calibration, every 15 minutes for observations in the very- 
extended configuration, and every 25 minutes for the compact 
configuration. We derived a flux density for J0423-013 of 3.3 Jy 
in July and 4.5 Jy in October and for J0334-401 the derived 
flux density was 1.2 Jy in July and 1.4 Jy in October. The 
quasar J0340-213, closer to NGC 1377 (1° away) but weaker 
(-0.4 Jy) than J0423-013 and J0334-401 was also observed dur- 
ing the very-extended configuration track for calibration- testing 
purpose: after the gain solution derived from J0423-013 and 
J03 34-401 was applied to J0340-213, it was found with no mea- 
surable offsets relative to its known position, and with a consis- 
tent flux, making us confident with the calibration of our very- 
extended configuration data. 

After calibration within the dedicated MIR/IDL SMA reduc- 
tion package, both visibility sets were converted into FITS for- 
mat, and imported in the GILD AS/MAPPING and AIPS pack- 
ages for further imaging. 

2.1. 12 CO 2-1 emission 

For the 12 CO 2-1 data, sets of visibilities were combined and de- 
convolved using the Clark method with uniform weighting. This 
results in a synthesized beam size of 0."65 x 0."52 with position 
angle PA=8°. We smoothed the data to a velocity resolution of 5 
km s -1 , yielding a Icr rms noise level of 16 mJy beam" 1 . 

Furthermore, data taken with the compact array were also 
deconvolved with natural weighting resulting in a synthesized 
beam size of 4."88 x 2."93 (PA= 176°)(lcr rms noise level of 
23 mJy beam -1 ). This provides a low resolution map to which 
we compare the high resolution data. 

The integrated intensity maps, velocity fields and dispersion 
maps are presented in Fig. [T] and the high resolution channel 
maps are shown in Fig. [2] 

2.1 .1 . The integrated 12 CO 2-1 intensity and molecular mass 

In the high resolution map (Fig.[T], top left) the molecular emis- 
sion is resolved and arises mostly from a region of 1" (100 pc 
at D=2l Mpc) in radius. (From a Gaussian fit we find a full 
width half maximum (FWHM) major axis of l."76 and minor 
axis l."4). Fainter emission extends out to a radius of 2-3". The 
integrated emission structure is quite complex with various fila- 
mentary features extending out from the center. A tilted X-like 
shape appears to emerge from the center. In the low resolution 
map the emission is mostly unresolved, but some faint emission 
is extending to the west and north-east from the central conden- 
sation. 

The integrated intensities and estimated molecular masses 
are presented in Tab.Q] The H2 column density towards the inner 
£>=58 pc (average beam diameter) is estimated to N(R2)= 3.4 x 
10 23 cm leading to an average nuclear gas surface density of 
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Fig. 1. Top panels: High resolution moment maps of 12 CO 2-1. Left: Integrated intensity with contour levels 2 x 
(1,2,3,4,5,6,7,8,9) Jy beam -1 kms" 1 and gray scale ranging from 2 to 20 Jy beam -1 kms" 1 . The cross marks the position 
(2:03:36:39.073 £:-20:54:07.08, which is the position of peak 7b( 12 CO 2-1). Centre: Velocity field with contour levels starting 
at 1713 kms -1 and then increasing by steps of 8.6 kms -1 . The grayscale ranges from 1713 to 1764 kms -1 . Right: Dispersion 
map where contours start at 5 kms -1 and then increase by steps of 5 kms -1 . Grayscale ranges from to 30 kms -1 and peak 
dispersion is 90 kms -1 . Lower panels: Low resolution moment maps of 12 CO 2-1. Left: Integrated intensity with contour levels 
1.3 x (1, 2, 4, 8, 16, 32, 64) Jy beam -1 km s" 1 and gray scale ranging from to 129 Jy beam -1 km s" 1 . Centre: Velocity field with 
contour levels starting at 1713 km s" 1 and then increasing by steps of 8.6 km s" 1 . The grayscale ranges from 1713 to 1764 km s" 1 . 
Right: Dispersion map where contours start at 6.4 kms -1 and then increase by steps of 6.4 kms -1 . Grayscale ranges from to 
40 km s -1 and peak dispersion is 65 km s -1 . Note that the high resolution moment 1 and 2 maps are on a somewhat different scale 
compared to the moment map to better show the details of the velocity field and dispersion map. 



S « 5 x 10 3 Mope" 2 . The average column density for the central 
390 pc is « 5 x 10 22 cm" 2 . Please see footnote to Tab.HJfor 12 CO 
to H2 mass conversion. Potential errors are further discussed in 
Sect. [3T31 

We investigated what fraction of the single dish flux was 
recovered with our interferometric observations by comparing 
with our IRAM 30m 12 CO 2-1 observations. We find a 12 CO 2- 
1 luminosity in the 9" IRAM beam of 3.5 x 10 7 K kms" 1 pc 2 
- to be compared to the SMA total 12 CO luminosity of 3.3 x 
10 7 K km s -1 pc 2 . Within the errors we conclude that all of the 
IRAM 30m flux has been recovered with the interferometer. 



2.1.2. Kinematics 

Velocity field: The high resolution velocity field (Fig. [T] top 
centre) shows complex dynamics dominated by red-shifted gas 
to the east of the center and blue- shifted gas to the west. Fitted 
velocity centroids range from 1690 to 1770 kms -1 . In the low 
resolution map, the position angle of the smoother velocity 
field is dominated by a PA^40° component and velocity cen- 
troids range from 1712 to 1764 kms -1 . The velocity centroids 
were determined through a flux weighted first moment of the 



spectrum of each pixel therefore assigning one velocity to a 
potentially complex spectral structure. 

Dispersion maps: The high resolution dispersion map (Fig. [T] 
top right) shows a distinct X-like shape - where the "x" itself 
has a PA of 35° - 40°. Dispersion (one dimensional, cr ld ) in 
the X-structure is typically 25 kms -1 , but - reaches 44 kms -1 
or higher in isolated regions. The large dispersion is caused 
by multiple, narrow, spectral features and/or broad wing-like 
components. The low resolution dispersion map shows a less 
complicated pattern where the dispersion is dominated by a 
bipolar large scale feature with peaks of <xid=40 kms -1 . The 
dispersion was determined through a flux weighted second 
moment of the spectrum of each pixel. This corresponds to the 
one-dimensional velocity dispersion (i.e. the FWHM line width 
of the spectrum divided by 2.35 for a Gaussian line profile). 

Spectra and channel map: The low resolution integrated 
12 CO 2-1 spectrum is characterized by narrow emission at 
the line center with FWHM line width £V^60 kms -1 and by 
triangular- shaped wings (or plateau) ranging from 1620 kms -1 
to 1860 kms -1 (see Fig. [3] top panel). Wings are therefore 
stretching out to ±120 kms -1 (projected). This line shape 
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Table 1. 12 CO 2-1 flux densities and molecular masses 

1.5 



Position^ (J2000) a: 03:36:39.076 (+ 0."2) 



Peak flux density c 

F\ux d 
(central beam) 
(whole map) 



6: -20:54:07.01 (+ 0."1) 
270+ ^(mJybeam" 1 ) 

20 + 0.05 (Jykms -1 ) 
169 + 0.5 (Jykms -1 ) 



Molecular mass e 
(central beam) 
(whole map) 



1.7 xlO 7 M 
1.5 xlO 8 M Q 



a) The high resolution 12 CO 2-1 data. Listed errors are Icr rms. 

b) The position of the 12 CO 2-1 integrated intensity. The peak T B is at 
0-03:36:39.073 £:-20:54:07.08 at V c =1734 kms -1 . 

c) The Jy to K conversion in the 0."65 x 0.52" beam is 1 K=14.5 mJy. 

d) A two-dimensional Gaussian was fitted to the integrated intensity 
image. 

e) The 12 CO luminosity L( 12 CO)=18.5 x D 2 x (<9) 2 x 7( 12 CO) 
K km s _1 pc 2 (where D is in Mpc, 6 in "and 7( 12 CO) in K km s _1 ) and 
M(H 2 )=3.47L( 12 CO) for Af(H 2 )//( 12 CO)=2.5 x 10 20 cm" 2 ). The conver- 
sion factor has been calibrated for 12 CO 1-0 emission- thus, if the 12 CO 
is subthermally excited the H 2 mass estimated from the 2-1 line should 
be corrected upwards. From our IRAM 12 CO 2-1 and 1-0 data we esti- 
mate that the 2-1/1-0 line ratio to 0.7 - and hence we correct out masses 
upwards by 30%. 



is also seen i n the SEST single-dish 12 CO 1-0 spectrum by 
iRoussel et al.l (120031) . The high resolution channel map reveals 
significant structure and complexity in the 12 CO 2-1 line 
emission per channel. Note that there is no evidence of a bright, 
broad compact emission from an unresolved nuclear source. 
Brightest line emission in terms of peak 7b (19 K) occurs at 
V=1734kms _1 which we adopt as the center velocity. 

Line wings: To identify the spatial origin of the line wings, 
we integrated the emission corresponding to velocities 1660 - 
1676 km s" 1 (blue wings) and 1813 - 1829 km s" 1 (red wings) 
(y sys =1734 km s" 1 ) in the low resolution 12 CO 2-1 data. The re- 
sulting map is presented as the top left panel of Fig. HI The blue 
and red wings are spatially separated by ^2" and with a PA of 
45°. In the lower left panel of Fig. |4] we present an image from 
the high resolution data. Here we have integrated the emission 
corresponding to the central channels: blue: 1718-1729 kms" 1 
and red: 1739-1750 km s -1 . In this figure, we see a velocity shift 
along an axis perpendicular to the one seen for the line wings of 
the low resolution data. 

The relative positional uncertainty is on the order of (beam 
size)/(s/n ratio). For the low resolution wing map this is less than 
0."7 and for the high res map it is 0."05 - 0."1 . Therefore the 
two orthogonal velocity gradients are observationally robust. 



2.2. 13 CO and C 18 2-1 emission 

For 13 CO and C 18 2-1, as the signal-to-noise ratio is lower, 
we used only data from the compact configuration, deconvolved 
using natural weighting and smoothed to a velocity resolution 
of 11 km s" 1 . The resulting beam size is 4."88 x 2."93 (RA. 
= 176°) and the Icr rms noise level is 13 mJy beam -1 . The 
low resolution 13 CO and C 18 spectra are presented in Fig. [3] 
The 13 CO is clearly detected and a factor of 9-10 fainter than 
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Fig. 3. Spectra of Top panel: 12 CO 2-1, Middle panel: 13 CO 2- 
1 and Lower panel: C 18 2-1. All three spectra are taken at po- 
sition a=03:36:39.07 £-20:54:07.00 in the low resolution data. 
13 CO is clearly detected and a factor of 9-10 fainter than 12 CO. 
C 18 is not detected and the 12 CO/C 18 2-1 peak line intensity 
ratio is >50. The fits to the 12 CO and 13 CO spectra are discussed 
in Sect. [3X3] 



12 CO. C 18 is not detected and the 12 CO/C 18 2-1 peak line 
intensity ratio is >50. The integrated intensity and velocity field 
maps are presented in Fig. [5] The 13 CO peak integrated emis- 
sion is offset by 0."9 (90 pc) to the south-east from the peak 
integrated 12 CO emission. Furthermore, the 13 CO velocity field 
appears more consistent with that of the high-resolution 12 CO 
data. 

2.3. Continuum emission 

We detect (4cr) weak 225 GHz continuum in NGC 1377 of 2.4 
(±0.45) mJy beam -1 in the compact array at the position of the 
12 CO 2-1 line emission. The source is unresolved in the 4." 88 x 
2."93 beam. 



3. Discussion 

3.1. A molecular outflow? 

We propose that the molecular emission of NGC 1377 is dom- 
inated by a disk-outflow morphology. The notion of a molecu- 
lar outflow is supported by the presence of 11 CO 2-1 line wings 
(Fig.\3§ and the spatial orientation of the wing-emission with re- 
spect to the line-center emission (Fig. [?]). The X- shape of the in- 
tegrated intensity and dispersion maps (Figs. \T\and\7^ is consis- 
tent w ith a biconical outflow structure. In addition jRoussel etaD 
(2003) find NIR H 2 (10) S(l) line emission towards the inner 
200 pc of NGC 1377 with an orientation of the velocity field 
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Fig. 2. High resolution channel map of the 12 CO 2-1 emission toward NGC 1377, as obtained from the combination of the compact 
and very-extended configuration data. Contours are drawn every 48 mJy beam -1 (3<x). The synthesized beam of 0."65 x 0."52 is 
shown in the bottom left corner. 



similar to the one we find for the 12 CO 2-1 line wings. Other 
alternatives to explain the kinematical structure of the 12 CO 2-1 
emission include a polar disk or inflowing gas. A polar disk is 
unlikely since this would require its velocity to be higher than 
the rotational velocity of the major axis disk. That the kinemat- 
ics is the signature of an outflow, instead of inflowing gas, is 
supported by the optica l dust absorpti on feature south of the nu- 
cleus visible in Fig. 1 in lRoussel et al] (120061) . The feature is spa- 
tially coincident with the blue- shifted line wing suggesting that 
the gas on the near side is coming towards us in an outflow. Note 
that non-circular motions in the galaxy plane may be present and 
will then add to the complexity of the dynamics. 



Below we will discuss the morphology and estimates of 
molecular mass and velocity of the suggested outflow, but we 
begin with a brief discussion of the nuclear disk. The fitted prop- 
erties of the molecular disk-outflow system are summarized in 
Tab. El 



3.1 .1 . Disk morphology and nuclear dynamical mass 

Disk morphology: To further investigate the notion of a disk- 
outflow system, we made moment maps of only the brightest 
emission (above 190 mJy beam -1 ) isolating the narrow line 
emission emerging from the inner rotating disk (Fig. |6). Note 
that the fitted PA of the disk of 110°-130° is consistent with it 
being perpendicular to the wing emission, the outflow, while 
the velocity field (right panel) is closer to a PA of ^95°. This is 
expected if an outflow is affecting spectral shapes and velocities 
even close to the nucleus of NGC 1377. The fitted disk diameter 
(Z)(FWHM)=0."60) should be viewed as a lower limit since 
we use only the top 23% of the emission. Disk inclination 
is roughly determined to be 45° - 70°. Uncertainties largely 
depend on the difficulty in small-scale separation between 
disk and outflow components. The stellar large scale disk is 
estimated to have an inclination of 60° (from the K-band image) 
and since it is consistent with the range of values for fits to the 
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Fig. 4. Left: Top panel: Blue (dark) (V= 1660 - 1676 km s" 1 ) 
and red (light) (V= 1813 - 1829 km s" 1 ) showing the location 
of the wings of the global 12 CO spectrum. This is done with the 
low resolution data for maximum signal to noise. Lower panel: 
Blue (dark)(V= 1718 - 1729 km s" 1 ) and red (light) (V= 1739- 
1750 km s" 1 ) part of the central part of the spectrum of the high 
resolution data. The velocities are selected to not overlap with 
the line wings. Contours are given every 3<x. Two velocity com- 
ponents - perpendicular to each other - can be discerned. Right: 
Top panel: Low resolution 12 CO 2-1 spectrum where the spec- 
tral regions in the line wings used to make the figure to the left 
are indicated. Lower panel: High resolution 12 CO 2-1 spectrum 
where the spectral regions in the line center used to make the 
figure to the left are indicated. 





Fig. 6. Integrated intensity and velocity maps of the high resolu- 
tion data where only emission in excess of 190 mJy was included 
to show the location and orientation of the nuclear disk. Position 
of cross is the same as in Fig.[T] Left panel: Integrated intensity 
where contour levels are 0.72 x (1, 3, 5, 7, 9) Jy beam -1 km s" 1 
and gray scale ranges from to 7.2 Jy beam -1 kms -1 . Right 
panel: Velocity field with contour levels starting at 1725 km s" 1 
and then increasing by steps of 5 kms -1 . The grayscale ranges 
from 1720 to 1750 kms . The synthesized beam is 0."65 x 
0. ,, 52. 



(where V rot is in km s" 1 and r in pc). This is only 2.3 times that of 
the estimated nuclear molecular mass implying a relatively high 
molecular gas mass fraction in the center of ^40%. A limit to the 
diameter of the molecular disk (which would be the extension of 
the nuclear disk) can be set to 200 pc based on the total extent of 
the SMA high resolution 12 CO 2-1 emission along the disk PA. 
This is of course sensitivity limited and a fainter, larger molec- 
ular extension to the disk is possible. The Keplerian dynamical 
mass of this 200 pc disk would be 1.3 x 10 8 M Q . 
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Fig. 5. Low resolution moment maps of 13 CO 2-1. Left panel: 
Integrated intensity with contour levels 0.84 x (1,3, 5, 7, 9, 11) 
Jy beam -1 kms -1 and gray scale ranging from 1 to 10 
Jy beam -1 kms -1 . Right panel: Velocity field with contour lev- 
els starting at 1717 kms -1 and then increasing by steps of 
8.6 km s" 1 . The grayscale ranges from 1713 to 1764 km s . 



nuclear disk we will adopt this value for the nuclear disk as well. 

Disk dynamics: The fitted maximum projected rotational ve- 
locity of the disk is 65 km s" 1 and for an inclination of 60° this 
results in a V rot of 75 km s~ . The dynamical mass of the 60 pc 
diameter nuclear disk is estimated to Md yn = 4 x 10 7 M Q using 
the Keplerian formula M dyn = 2.3 x 10 8 x (V rot /100) 2 x (r/100) 



3.1 .2. Outflow morphology 

The X- shaped morphology of the high resolution integrated in- 
tensity and dispersion maps (see Fig. [7]) suggests the molecular 
gas outlines a symmetric biconical structure centered on the nu- 
cleus. The cone has an opening angle of 60° - 70°. Such a large 
opening angle could explain the complex velocity pattern ob- 
served in the high resolution map. The peak intensity and peak 
dispersion are roughly correlated in space. The large dispersion 
may be caused by line-of- sight effects in maximum path-lengths 
along the cone walls, and also by multiple spectral components 
where extended disk emission overlaps with emission from the 
outflowing gas. 

Note, however, that the molecular outflow may not be 
well-collimated or in a simple well-ordered outflow pat- 
tern, but instead filamentary and patchy and the molecu- 
lar emission is likely a mix of clumpy and extended gas. 
(Properties of wind-cone struc ture s hav e been modeled by 
iTenorio-Tagle & Munoz-Tunonl (e.g. [l998l) ). The X- shape struc- 
ture also appears to be evident in the dust morphology of 
NGC 1377. The optical dust absorpti on feature south of the nu- 
cleus (Fig. 1 in iRoussel et al.l (120061) ) shows a V- shape which 
is spatially coincident with the blue part of the molecular out- 
flow. The red- shifted molecular outflow is then hidden behind 
the galaxy with the blue-shifted outflow projected in front of the 
disk on the southern side. 
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3.1 .3. 12 CO luminosity and molecular mass of the outflow 

Without a detailed model the 12 CO 2-1 luminosity in the out- 
flow is difficult to determine. Below we adopt two approaches to 
obtain upper and lower limits to the luminosity. 
When estimating a lower limit to the luminosity in the outflow 
we restrict ourselves to the (1628 - 1655 kms" 1 ) and (1813 
-1840 kms -1 ) part of the plateau component (i.e. 79 kms -1 
~|V-y sys |n06kms- 1 )to ensure that the rotational components 
are not included. We find a flux of 13 ± 0.3 Jy beam -1 . Since we 
are only considering a narrow (27 km s -1 ) part of the 120 km s -1 
outflow this will be an underestimate of the total 12 CO 2-1 flux 
in the outflow. If the opening angle is as large as we discuss in 
Sect. 13.1.21 there will be significant projected outflow emission 
also at systemic velocities. 

To obtain an upper limit of the outflow luminosity we fit a 
Gaussian to the line wings in the low resolution spectrum (see 
Fig©. We first fit a Gaussian to the 13 CO 2-1 spectrum (as- 
suming that the 13 CO has only a small contribution from the 
wings) where all parameters are free. We obtain: S pe ak=0. 15+0.1 
Jy beam" 1 ; V c =1740+2 kms -1 and £V(FWHM)=61±5 kms -1 . 
Then, for the 12 CO 2-1, we fit two Gaussian curves, using the 
parameters from the 13 CO fitting. The narrow component has 
Speak= 1.6+0.05 Jy beam -1 with V c and SV fixed. The broad 
component in the wings has, with V c fixed, 5 , p e a k=0.4+0.05 
Jy beam" 1 and £V(FWHM)=151+8 kms" 1 (and the full width 
of the line wings is +120 kms -1 (projected)). The division of 
flux between the line core and outflow component is 104 and 
65 Jy beam -1 km s -1 respectively. Since we detect 20 Jy beam -1 
kms -1 in the central high resolution beam, the narrow compo- 
nent flux is mostly located outside of the inner 0."6. This flux 
may be distributed as an extension to the nuclear disk we dis- 
cuss above in Sec. 13.1.11 

To estimate the molecular mass in the outflow we adopt a stan- 
dard 12 CO-luminosty to M(H2) conversion factor (see footnote 
to Tab. Q]). This results in M(H 2 ) ut=l.l * 10 7 M Q for the lower 
limit and for the upper limit M(H 2 ) ut=5.4 x 10 7 M Q . The valid- 
ity of a standard 12 C O to M(H2) conversion factor in galaxies is 
heavil y debated (e.g. Narayanan et al.l l2012: Wad a & Tomisakal 
120051: IPagfione et al.1 l200lh . The global 12 CO/ 13 CO 2-1 line 
intensity ratio of 9-10 suggests that the molecular gas in 
NGC 1377 is not optically thin (but 13 CO limits in the line wings 
are not good enough for a reliable line ratio for the highest veloc- 
ities in the outflow). It is also unclear whether the interferometer 
would pick up emission that is optically thin since it tends to 
have low surface brightness. We therefore suggest that applying 
the conversion factor here does not overestimate gas masses by 
more than factors of a few. 

3.1 .4. Mass outflow rate, velocities and acceleration. 

The projected peak outflow velocity is -120 kms -1 (sec- 
tion !2.1.2l) and if the opening angle of the outflowing gas is -60° 
(and the disk inclination is /=60°) then the actual maximum out- 
flow velocity is -140 kms -1 . The outflow has reached out to 
-200 pc (corrected for opening angle and inclination) which, for 
y out =140 km s" 1 gives t dyn = R/V=1A x 10 6 yrs. The molecular 
mass in the outflow is M out = 1.1 -5.4x 10 7 M Q (see Sect. 13.1.31) . 
Thus, the molecular mass loss rate SM/St = M- 8-38 M yr" 1 . 

We do not know the escape velocity of NGC 1377 hence 
it is not possible to know whether the gas will leave t he sys- 
tem as a whole, or just the nuclear region of NGC 1377. Martin 



Table 2. Summary of outflow and disk properties^. 



Outflow 




A/Tolppnlar mflss' 


1 1 - 5 4 x 1 7 


Extent: 


200 pc 


Orientation: 


PA=35° - 45° 


Opening angle: 


60° - 70° 


Age: 


1.4 x 10 6 yrs. 


Outflow velocity: 


140 kms" 1 


M: 


8-38M yr 1 . 


Energy: 


2 - 10 x 10 54 ergs 


Power: 


0.4 - 2 x 10 41 ergs s 


Nuclear disk 




Diameter: 


0."60 (60 pc) 


Orientation: 


PA=110°-130° 


Rotational velocity: 


75 kms" 1 


Dynamical mass: 


4 x 10 7 M Q 


Molecular mass: 


1.7 xlO 7 M 


H2 column density: 


3.4 x 10 23 cm" 2 


H 2 surface density: 


5 x 10 3 M Q pc~ 2 



a) Properties are derived and discussed in sections 12.1.21 ; 13. 1.T1 ; I3. 1.21 ; 
15X51 : 15X51 



(2005) sets a limit for the escape velocity as Vesc=3V ro t and if 
the rotational velocity is 75 km s -1 then the gas will not escape 
NGC 1377. 

An inspection of the central spectrum reveals that broad fea- 
tures are present in the very centre as well suggesting that the 
gas is accelerated close to the nucleus (see lower right panel 
of Fig. |4]). Compared to outflows in M 82 and NGC 253 it ap- 
pears that the acceleration to the final velocity happens on a short 
length scale - smaller than 30 pc. 

3.2. Energy and power of the outflow 

We use the prescription by Veille uxet al.l (1200 ll) to calculate the 
kinetic energy in the outflow: E^i n = £buik+ ^turbulent- This kinetic 
energy represents a lower limit to the energy. The work of lifting 
the gas out of the nucleus should be included but since we do not 
know the actual gravitational potential depth of the galaxy this 
term cannot be added at this stage. If we assume that the molec- 
ular gas in the outflow is moving at constant speed y out =140 
km s" 1 we obtain E buXk =i: i (l/2)Sm i Vi 2 =(l/2)M 0Ut x V 2 ut ~2x 10 54 
ergs - when adopting the lower limit to M out . We can however 
only measure an upper limit to the turbulence of the gas of 
<Tid=25 kms -1 (due to effects of spatial resolution, rotation and 
the outflow) and since this term is already smaller than the con- 
tribution from the bulk motion we do not include it and in total 
we estimate E^ n to ~2 x 10 54 ergs. For the upper limit to M out 
we get i^ot^iO 55 ergs. The power in the outflow is 0.4 - 2 x 10 41 
ergs s" 1 which corresponds to a luminosity of 1 -5 x 10 7 L©. The 
ratio between the wind luminosity and the FIR luminosity is then 
2 x 10" 3 - 10" 2 . From now on we adopt the more conservative 
lower value of the luminosity and mass outflow rate. 

3.3. What is driving the outflow? 

The unusual and obscured nature of NGC 1377 makes it difficult 
to find an obvious source of the molecular outflow. Is it a buried 
AGN or a young starburst? Below we investigate some possible 
scenarios. 
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3.3.1 . Winds from supernovae and massive stars 



The limit on L(1.4 GHz) is 5.3 x 10 19 W Hz" 1 dRoussel et al.1 
2003| ). This corresp onds to an upper limit to the SNR of 5 x 10" 4 



yr 
34 in 



Condon 1992). Inserting this limit into equations 10 and 
Murray et aTT(l20Q5h we find that SNR fall short by an or- 
der of magnitude to drive the outflow - by momentum or energy. 
The lack of radio continuum from SNR is consistent with th e 
nondetection of NIR [Fe II] emission by iRoussel et al.l d2Q06h . 
However, an outflow may also be driven by the wind-momentum 
of massive stars. The total momentum in th e NGC 1377 outflow 
is 1 .5 x 10 9 M Q km s" 1 . Using Starburst99 dLeitherer et al.ll 19991) 
we find that a 10 6 M cluster with SFR of 1/3 M Q yr" 1 has a wind 
momentum of 10 7 M km s -1 when integrated over 3 Myrs. Thus 
an SFR of ~30 M yr -1 would be required. From the FIR lumi- 
nosity we estimate a SFR of 0.4 M© yr" 1 for NGC 1377 (assum- 
ing th at all of the FIR is due to star formation) and Roussel et al.l 
( 2003) estimate a rate of 1.8 M© yr -1 . Thus it appears that the 
SFR is too low by about one order of magnitude to push out the 
molecular gas by the wind momentum of massive stars. 

Could the outflow be driven by internal ram pressure from 
hot gas? Only upper limits to the thermal free-free emission in 
the radio has been obtained. In Fig. 8 of IRoussel et al.l (120031) 
predictions of the radio synchrotron and free-free emission based 
on the IR are plotted. For NGC 1377 the VLA 1.4 GHz upper 
limit is about a factor of 4-5 below the expected value. There is 
thus no direct evidence of much hot ionized gas. Both the Ha 
and [N II] optical emission line s are w eak and Bry and Paa are 
undetected (IRoussel et al.ll2003Ll2006h . This may be an effect of 
the high extinction so to complete the picture of the hot gas of 
NGC 1377, X-ray observations will be essential. 

Note also that the symmetry of the wind is inconsistent with 
the idea of external ram pressure driving the wind, as is the pres- 
ence of NGC 1377 in a low density galaxy group where ram 
pressure forces are expected to be low. 



3.3.2. Radiation pressure 



Mur ray etall (120051) suggest that radiation pressure (from both 
star formation and AGNs) from the continuum absorption and 
scattering of photons on dust grains may be an efficient mech- 
anism for driving cold, dusty gas out of a galaxy. The lack 
of evidence for supernovae and the presence of a dusty, com- 
pact IR source make s this an attractive scenario for NGC 1377. 
Mur ray etall (120051) suggest that for a radiative pressure driven 
outflow the momentum flux MV should be comparable to that in 
the radiation field L/c. For NGC 1377 the outflow momentum 
flux exceeds L/c by factors of a few which is within the uncer- 
tainties for outflow structure, velocity and mass. We suggest that, 
with current information, radiation pressure is a potential driv- 
ing mechanism for the outflowing gas in NGC 1377, but the fit 
is far from perfect and other processes cannot be excluded. 

3.3.3. AGN or starburst? 



We use the iGraham et all (l201lb (their Figs. 2 - 4 and Tab. 2) 
calibration of the stellar velocity dispersion and BH mass to 
estimate th e mass of a superm assive black hole (SMBH) in 
NGC 1377. IWegner et all (120031) list a value of the stellar ve- 
locity dispersion of <x=83 kms" 1 for NGC 1377 in their on- 
line database. This gives a SMBH mass in the range of about 
4 x 10 5 - 6 x 10 6 M Q with a most likely value of about 1.5 x 10 6 
M©. If all the FIR emission observed would be due to accretion 
onto the SMBH it would then operate at -10% of its Eddington 




Fig. 7. Overlay of high resolution grayscale dispersion map 
(grayscale ranging from to 34 kms -1 ) on contours of 
integrated intensity with contours 2 x (1,2,3,4,5,6,7,8,9) 
Jy beam -1 kms -1 . Note how regions of high dispersion form 
an X- shaped structure generally aligned with a similar structure 
in the integrated intensity map. 



limit. A small acceleration region is consistent with a compact 
source of energy and momentum. Radiation pressure accelera- 
tion minimizes hard shocks which results in low radio emission 
and acts on dense gas, consistent with the IR H2 emission and 
the suggested 12 CO 2-1 molecular outflow. It is also easier to 
bury a compact energy source behind the gas and dust, giving 
rise to the deep silicate absorption. Strong FIR emission lines 
from ionized gas in an AGN are expected to be highly localized 
around the SMBH, and thus are easier to hide. 

In contrast, a dust-enshrouded nuclear starburst needs to be 
located well inside a radius of 29 pc (to launch the outflow). The 
current molecular mass here is estimated to 1 .7 x 10 7 M Q , 40% 
of the dynamical mass. IRoussel et al] (120061) suggest a starburst 
with a mass 1.6 x 10 7 M assuming an age of 1 Myr, a Salpeter 
IMF between 0.1 and 120 M Q and 3 x 10 4 O-stars to produce 
the bolometric luminosity of NGC 1377. The star formation ef- 
ficiency (SFE) of the starburst then had to be at least 50% and 
all of its Paa, and Bry emission has to be absorbed. This implies 
that the stars must be deeply buried inside the nuclear dust cloud 
putting even stricter limits on their radial distribution. In addi- 
tion, the starburst scenario requires that almost all of the stellar 
mass in the center of NGC 1377 has to be very young. 

3.4. Comparing with other molecular outflows 

Outflowing molecular gas has been detected towards several 
galaxies of various luminosities. In some cases the existence of 
an outflow has been infe rred from 12 CO line wings - such as 
for the LIRG NGC 3256 dSakamoto et al.ll2006h . In some cases, 
the molecular outflow was discovered in al ready iconic opti- 
cal outflows, such as is the case for M 82 (Na kai etall 1987: 
Walter et al. 2002). There is also mounting evidence of massive, 
V> 1000 kms" 1 outflows in AGN/starburst driven ULIRGs (e.g . 
Sak amoto et all 120091: iFeruglio et al.l l2010t IChung et al.l l2011h 
Stur m et aljl201lH Aalto et al.ll2012l) . 



S. Aalto et al.: CO in NGC 1377 



9 



NGC 1377 can also be compared to the AGN-driven molec - 
ular outflow of the SO galaxy NGC 1266 (lAlatalo et al.ll2011l) . 
Nuclear gas depletion time scales are similar to those estimated 
for NGC 1377, and both galaxies have compact, r < 100 pc nu- 
clear disks confining the outflowing wind. NGC 1266 is more 
FIR luminous than NGC 1377 with an order of magnitude 
more molecular gas and its molecular disk has a factor of five 
greater gas surface density. Also, NGC 1266 is not radio defi- 
cient (g=2.3) unlike NGC 13 77. Estimating the p ower of a radio 
jet from the total radio flux, lAlatalo et aiT(l201 lb attributed the 
molecular outflow to the mechanical w ork of the radio jet of the 
AGN (see also Matsush ita et al.l {2007) for a similar case). 

What sets NGC 1377 apart is its unusual IR properties and 
extreme q suggesting that the nuclear activity is in a transient 
phase of its evolution or that NGC 1377 is exceptional in some 
other aspect. This is illustrated by the fact that observationally 
the q parameter is in a strikingly small range in large statistical 
samples of galaxies. In particular, the upper-bound of q is very 
clear and a > 3 is really exceptional (see for example Fig. 6 of 
lYun et al] d200lb . 

Comparing the energy and mechanical luminosity of the 
NGC 1377 outflow to one AGN-driven (NGC 1266) and one 
starburst driven (NGC 3256) outflow we find them to be similar 
(within factor of a few). The ratio between the wind luminosity 
and the FIR luminosity Lmech/^FiR, however, for NGC 1377 is 
similar to, or higher, than that of NGC 1266, and at least one 
order of magnitude higher than that of the starburst driven out- 
flow of NGC 3256. More studies will show if this difference is 
related to the underlying driving mechanism. The L m ech is calcu- 
lated only on the molecular flows in all three galaxies. 

3.5. Nuclear gas properties 

3.5.1 . Is there a Compton-thick nuclear dust and gas cloud? 

For NGC 1377 to hide its power source also at X-ray wave- 
lengths large columns of absorbing material has to reside in 
front of the nucleus. We estimate the average A/^Fb) in the in- 
ner r=29 pc of NGC 1377 to be 3.2 x 10 23 cm" 2 . (This corre- 
sponds to an Ay of 290 (for N n (cm" 2 ) = (2.21 ± 0.09) x 10 21 A V 
(iGiiver & O zel 2009).) For an absorber to be Compton-thick an 
A^h>10 24 cm -2 (Ay>450) is required. Thus, a Compton-thick ab- 
sorber, if any, should therefore be smaller than our beam. If the 
12 CO luminosity of our 12 CO 2-1 beam is located in a region 
with a radius 1 17 pc the resulting Nu would be ! 10 24 cm -2 . 



3.5.2. Gas physical conditions 

We find a peak 12 CO 2-1 temperature of 19 K in a 58 pc beam 
which serves as a strict lower limit to the gas kinetic temperature 
in this region. The dust temperature is difficult to determine but a 
fit to the 25-100 jim data result s in a black body te mperature of 
80 K with a diameter of 37 pc dRoussel et alj|2003h . Comparing 
our 1 mm continuum flux to the FIR SED suggests it is consistent 
with a v 3 slope - either implying fi= 1 - or that the dust is already 
opaque at A=100 //m. This possibility, if confirmed, would be 
consistent with the dusty opaque core suggested from the deep 
silicate absorption. There is no evidence of an excess of 1 mm 
flux emerging from cold (10-20 K) dust. 

We can consider three simple scenarios for the physical condi- 
tions of the nuclear molecular gas: 

1. The molecular gas is cold and optically thick. The 12 CO 2- 
1 emission is optically thick, fills the central beam and its 



brightness temperature traces the kinetic temperature of the 
gas at 29 pc radius from the nucleus. This requires a very 
steep, radial temperature gradient to allow for the buried 
warm dust. 

2. The molecular gas is optically thin In this scenario, the 12 CO 
2-1 brightness temperature does not reflect the kinetic tem- 
perature of the gas even if it fills the central beam and will 
instead be oc T QX r. Given the large columns of absorbing gas 
and dust this scenario seems unlikely. 

3. The molecular gas is warm and optically thick A molecu- 
lar medium of warm gas clumps that only fills a fraction of 
the beam. Alternatively, the 12 CO peak is unresolved in the 
beam. If we assume that the emission is emerging from the 
same region as the inferred dust diameter (of 37 pc) the peak 
12 CO surface brightness temperature would be 50 K. 



The above scenarios assume a single-component ensemble 
of clouds. A more realistic view on the nuclear molecular ISM 
is that it likely consists of both dense in > 10 4 cm -3 ) gas clumps 
embedded in a lower density (/? = 10 2 - 10 3 cm " 3 ) molecular 
medium (see discussion in (e.g. lAalto et al.lll994l) ). The 12 CO 2- 
1 brightness will then reflect the relative filling factor of the low- 
and high density gas. The detection of bright HCN 1-0 emission 
(llmanishi et alJl2009b towards the inner 2"of NGC 1377 sug- 
gests gas densities >10 4 cm -3 (unless the HCN is IR pumped). 

Is the molecular gas concentr ation forming sta rs? From the 
Kennicutt-Schmidt (KS) relation dKennicutl fl998) the expected 
SFR from the nuclear gas concentration is 0. 1M© yr" 1 . This falls 
short by factors of 4 - 1 8 to explain the inferred SFR if all FIR 
emission is emerging from star formation. However, the KS re- 
lation at r<100 pc may should have a large scatter since star 
formation i s unlikely to be a ste ady- state process on such small 
scales. (e.g. lOnodera et al.l2010h . There should be some star for- 
mation going on even if it is difficult to infer the actual SFR. 



3.6. The evolutionary implications for NGC 1377 



3.6.1 . Origin of the molecular gas 



Just as is the case for NGC 1266 (Alatalo et al. 2011), one ques- 
tion for NGC 1377 is how the gas got to the very center of the 
galaxy so efficiently. So far, no clear evidence of an interaction 
has been presented. The range of possible PAs of the nuclear 
disk includes a difference by about 40° compared to the stellar 
disk of NGC 1377. An inclination difference between the gas 
and stellar disks wou ld indeed suggest an external source for the 
gas ( iDavis et al]|201ll) . It is interesting that the evidence so far 
suggests that most of the gas has ended up in molecular form. 
If NGC 1377 accreted a late type dwarf galaxy or gas from the 
outer regions of a large disk galaxy, the gas is more likely to be 
atomic. Somehow the gas must have become shocked and com- 
pressed along the way. However, the lack of 1 .4 GHz radio con- 
tinuum emission makes it difficult to assess the amount of cold 
atomic hydrogen in the center of NGC 1377 - which normally 
could be tested by searching for absorption in the 21 cm HI line. 

The gas concentration may instead be the result of slow, sec- 
ular evolution such as accretion along a bar (e.g. I Jo gee et all 
l2005h . very slow exte rnal accretion (galaxy harassment), or stel- 
lar return (e.g. IWelch & Sage 2003). All of these - or combina- 
tions of them - are possible scenarios, but the problem remains 
of getting such a concentration of molecular gas to the core of 
NGC 1377. 
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3.6.2. The nuclear activity 

Observing the molecular properties of NGC 1377 allow us to 
probe the early stages of nuclear activity and feedback mecha- 
nisms in active galaxies. Depending on what is driving the FIR 
emission and outflow of NGC 1377 we are either witnessing the 
growth and feeding of a black hole in an early type galaxy - a 
process fueled by the late influx of molecular gas. Alternatively, 
the gas fuels the early triggering phase of a nuclear starburst in- 
creasing the stellar mass instead of the central black hole. 

Another important key to the evolutionary stage and nature 
of the activity is the age of the outflow. If the SMA data recov- 
ers the full extent of the outflow, then this sets an upper limit to 
its age of 1.4 Myrs. High sensitivity observations of 12 CO 1-0 
as well as deep imaging of dust absorption features, Ha and a 
search for HI will help determine if the outflow extends beyond 
the 12 CO 2-1 emission. Such evidence would further constrain 
the age of a buried power source. An implicit assumption here 
is that this activity is a one-time event (i.e., not recurrent in a 
short time scale). The ongoing, mass outflow rate would clear 
the central region of gas within the short time period of 5 - 25 
Myrs. The growth of the SMBH and/or nuclear star formation 
will cease within this time period. However, if the gas is unable 
to escape the galaxy we might be observing a failed wind (e.g. 
Dorodni tsyn et al.ll201 ll) . The gas may then fall back onto the 
galaxy again where the returned gas serves as a reservoir for 
future activity. Further studies will address this possible evolu- 
tionary scenario. 

4. Conclusions 

We have imaged 12 CO and 13 CO 2-1 in the FIR-excess galaxy 
NGC 1377 with the SMA: 

1. We find bright, complex 12 CO 2-1 line emission from the 
nuclear region of the FIR-excess SO galaxy NGC 1377 the 
structure of which is reminiscent of a disk-outflow system. 
The 12 CO 2-1 line wings, the spatial orientation of the 
wing-emission with respect to the line-center emission, 
its correspondence with the optical dust features and the 
X-shape of the integrated intensity and dispersion maps 
support this interpretation. Estimated outflow parameters 
are: M 0Ut (H 2 )>l x 10 7 M ; M>8 M Q yr" 1 ; V ou t>140 kms" 1 ; 
extent^200 pc; opening angle=60° - 70°; PA^40° and an 
age of ~1.4xl0 6 yrs. 

2. We suggest that the age of the outflow supports the notion 
that the nuclear activity is young - a few Myrs. The dusty 
center of NGC 1377, the lack of significant numbers of re- 
cent supernovae, or of a dominant hot wind, may imply that 
the outflow is driven by radiation pressure from a compact, 
dust enshrouded nuclear source. However, other driving 
mechanisms are possible. Accretion onto a 10 6 M SMBH - 
i.e. a buried AGN - is a possible power source. Alternatively, 
the photons may emerge from a very young (1 Myr), 
pre-supernova starburst, but we find this explanation less 
likely due to the extreme requirements of compactness and 
youth of the source. The implied mass outflow rate is large 
enough that it will clear the nuclear region of gas within 
5-25 Myrs and will, at least temporarily, shut off the nuclear 
activity and growth. If the outflowing gas rains back onto 
the galaxy, however, it may serve as fuel for a future period 
of nuclear activity and/or star formation. 



3. Molecular masses are estimated through the adoption of 
a standard 12 CO to H 2 conversion factor - with all the 
implied uncertainties. In the whole 12 CO 2-1 map we detect 
M(H 2 )=1.5 x 10 8 M and in the inner r=29 pc M(H 2 ) is 
estimated to 1.7 x 10 7 M with an H 2 column density of 
Af(H 2 )= 3.4 x 10 23 cm -2 . This is not large enough for the 
core to be Compton thick and a Compton-thick absorber, if 
any, should therefore be smaller than our beam. The aver- 
age gas surface density within r=29 pc is X=5 x 10 3 M Q pc -2 . 

4. The 12 CO 2-1 brightness temperature sets a lower limit to 
the nuclear gas kinetic temperature of 19 K, but if cloud fill- 
ing factors are less then unity then gas temperatures will be 
significantly higher. 
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